Spin-transfer torque on a single magnetic adatom 
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We theoretically show how the spin orientation of a single magnetic adatom can be controlled by 
spin polarized electrons in a scanning tunneling microscope configuration. The underlying physical 
mechanism is spin assisted inelastic tunneling. By changing the direction of the applied current, the 
orientation of the magnetic adatom can be completely reversed on a time scale that ranges from a 
few nanoseconds to microseconds, depending on bias and temperature. The changes in the adatom 
magnetization direction are, in turn, reflected in the tunneling conductance. 

PACS numbers: 



There is now a fast growing interest in controlling 
the spin orientation of a single or few magnetic atoms 
in a solid state environment for future spintronics and 
quantum information devices. So far this has only been 
achieved using optical methods When a current 

flows through a magnetic region it becomes spin polar- 
ized due to exchange coupling between transport elec- 
trons and localized magnetic moments. The back action 
of transport electrons on the magnetic moment, known as 
spin-transfer torque [2|, can be used to rotate the mag- 
netization in nanopillars made of billions of atoms Q. 
The magnetization of such a large number of atoms can 
be described with a single classical vector and the cur- 
rent driven magnetization switching is properly mod- 
eled by Landau-Lifsthitz equations extended with the 
spin transfer term proposed by SlonczewskiQ. Current 
induced magnetization switching has been reported in 
much smaller nanomagnets, made of 100 atoms but 
still in the semiclassical domain. 

In the present letter we theoretically show how a spin 
polarized current can be used to manipulate the spin 
state of a single transition metal atom deposited on a 
insulating monolayer on top of a metallic surface [see 
Fig. HJa)]. The proposed implementation combines two 
alternative strategies used so far to probe the spin of a 
single atom on a surface: inelastic electron tunnel spec- 
troscopy (IETS) with non-magnetic tips on one side and 
a spin polarized tunneling current on the other. The 
first technique has demonstrated that conveniently iso- 
lated Mn, Fe and Co adatoms have quantized spin an- 
gular momentum along a well defined magnetic easy 
axis [E, S 0, S] • These experiments also demonstrate that 
tra nsp ort electrons and local spin are exchanged coupled 

[ToL [Tl]|. In contrast, experiments with ferromagnetic 
tips are based upon the sensitivity of conductance to the 
relative spin orientation of the tip and adatom (magne- 
toresistance) (l2| . 

Here we present a fully quantum mechanical theoreti- 
cal analysis showing that, under the influence of spin po- 
larized tunneling current, the spin of a single magnetic 
atom can be directed either parallel or antiparallcl to the 
magnetic moment of the tip (or the surface). When the 



current is spin polarized along a given direction ft, be- 
cause cither the tip or the substrate are ferromagnetic, a 
fraction of the tunneling electrons exchange one unit of 
spin with the magnetic adatom. This induces a net trans- 
fer of spin along n towards the magnetic atom. This cur- 
rent driven spin torque competes with the adatom spin 
relaxation provided by its weak coupling to the tip and 
substrate. The sign of the transfer is determined by the 
direction of the current whereas the efficiency is deter- 
mined by the magnitude of the current. 
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FIG. 1: (Color online), (a) Scheme of the proposed setup: 
a magnetic STM tip and a magnetic adatom on a insulating 
monolayer deposited on a metal, (b) Two of the spin assisted 
tunneling events of Eq. ([2]). Spin transitions when electrons 
flow from tip to surface ((c) and upper diagram in (b)) and 
from surface to tip ((d) and lower diagram in (b)). Dashed 
transitions (III) and (IV) are less efficient when the tip is 
polarized (see text). 

In order to study the spin dynamics of the magnetic 
adatom under the influence of the spin polarized tunnel- 
ing electrons we derive and solve the master equation for 
the eigenstates \M) of the single ion spin Hamiltonian 



H s = DS' z 2 + E(S' x 2 -S' y 2 ). (1) 
Thus, the local spin is described quantum mechani- 
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cally, in contrast to previous works 13|, [l4| . This is cru- 
cial to account for the IETS experiments Q. Here the 
prime denotes that the spin quantization axis is chosen 
with z perpendicular to the surface. For convenience, 
we choose the quantization axis of the all the spins in 
the Hamiltonian along ft which makes it necessary to ro- 
tate 7^s when ft is not perpendicular to the surface. The 
value of the total spin S 2 , and the magnetic anisotropy 
coefficients D and E change from atom to atom and also 
depend on the substrate [1, 0] . 

The scattering rates between states with different spin 
projection \M) arise from the exchange interaction of the 
adatom spin to delocalized electrons in tip and surface, 
hereafter denoted as the reservoirs and labeled with the 
indeces rj = T, S. The total Hamiltonian of the adatom 
coupled to the reservoirs reads H — Htip+Hsm+'Hs + V . 
The first three terms describe the tip, surface and adatom 
as decoupled systems, whereas the V term introduces in- 
teractions between them. The first term describes the 
electrons in the ferromagnetic tip magnetized along the 
direction n, a static parameter in our theory. The second 
term is the Hamiltonian of the non-magnetic electrons of 
the metallic surface. Thus, the quantum numbers of the 
transport electrons are their momentum k, spin along the 
n axis, er, and reservoir 77 = tip, sur. With this notation 



we write Ji 



tip 



yn( k ) c kar, c k<yn • Since 



we consider a non-magnetic surface, we have f. a ,s{k) 
es{k). All the results of this paper are trivially gener- 
alized to the case of a non-magnetic tip and a magnetic 
surface. 

The coupling of the atomic spin S and the conducting 
reservoirs has the form [1 51 ] 



V= Tjjtf^aQi, k ) 7^ S a C k(Jrl Ck>cr'r)' > (2) 

a,fc,fc' ,<t.(t' .r\.r\' 



where the index a runs over 4 values, a — x,y,z, and 
a = 0. We use r^ a ' and S a for the Pauli matrices 
and the spin operators in the n frame, while Sq = I 
is the identity matrix. T VJ1 > _ a (k,k') for a — x,y,z is 
the exchange-tunneling interaction between the localized 
spin and the transport electrons and potential scattering 
for a — 0. Attending to the nature of the initial and final 
electrode, Eq. © describes four types of exchange inter- 
action, two of which contribute to the current, the other 
two conserving the charge difference between tip and sur- 
face. The former are crucial to account for the magnetic 
IETS Q and tend to "heat" the spin of the adatom. 
The other two provide an efficient cooling mechanism, 
through a Korringa like spin relaxation, and were not 
included in previous work 

BE, EI- Following Ander- 
son [yl, we assume that (|2|) arises from kinetic exchange. 
The momentum dependence of T^ jl; ' iQ (fc, k') can have im- 
portant consequences in the conductance profile [TtJ in 
a energy scale of eV, so it can be safely neglected in 
IETS. We thus parametrize T n , n ', a =o(k, k') — v^v^Tq 



and TJj^/ (fc, k') = v^v^Ts, where Tg is the same for 
a = x,y,z and u sur and Vf lp are dimensionless factors 
that scale as the surface-adatom and tip-adatom hopping 
integrals. This parametrization sets a relation between 
spin torque and spin relaxation. 

The occupation of the spin states \M), Pm, are gov- 
erned by the master equation 
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M,M' ' (3) 



M ' .r\r\' 



where W^m' are the scattering rates from state M to 
M' induced by interaction with a quasiparticle which is 
initially in reservoir rj and ends up in r/. Eq. ([3]) docs 
not include spin coherences. This approximation is good 
provided that spin decoherence is faster than spin relax- 
ation, which is known to be the case due to hyperfine 
coupling [l[ in Mn atom. After some algebra, the scat- 
tering rates can be written as: 



W n ^ n - 

M,M' " 



-G(Am,M + V n - AV) S M M" ( 4 ) 



where Q{oS) = u> (1 — e~> ) are the phase space fac- 
tors associated to quasiparticle scattering, Am,m' = 
Em — Eftjr is the energy change of the atom, fj> v is the 

chemical potential of electrode rj and YP^J 1 M , are spin ma- 
trix elements 
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where Sf' M ' = {M\S a \M') and p 7ha is the density of 
states at the Fermi energy for spin a in the electrode rj. 
These equations show that, for a ferromagnetic tip, spin 
flip and spin conserving rates are different, in contrast to 
the case of a non-magnetic tip [9( . The rates that induce 
changes in the adatom spin population are all propor- 
tional to To, and can be classified in three groups: intra- 
tip, intra-surface and tip-surface rates (oc v^ ip , tx i£ ur , 
and oc v+^v 2 ^, respectively). In the first two, the capa- 



^tip^sur 

bility to transfer energy to the magnetic atom is given 
by the temperature, being the release of energy from the 
atom to the electrodes always allowed. In contrast, the 
rates carrying current can transfer energy to the atom 
even at zero temperature, provided that there is a bias 
voltage eV — /ig — fir- 

On top of the T| rates, there are scattering processes 
of order Tq and TqTs that do not change the occupa- 
tions but contribute to the current. Evaluated to the 
same second order in V than the rates in the master 
equation, the expression of the current has three terms, 
I = Iq + I air + 1 in- The first term is elastic and spin 
independent. The second is the elastic but sensitive to 
the relative spin orientation of the adatom and the tip 
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magnetic moment. The third term comes from the in- 
elastic exchange processes being proportional to Tj. The 
expressions for the three contributions to the current, 
analogous to those derived elsewhere [Hj|. are: 

Io + Imr = -~G a (1 + x{S z )V T ) i-(-eV)) (6) 

*™ = — £ [i-(*M,M>-eV)J2\^' M '\ 2 

6 M,W a 

+V T i+(A M ,M> - eV)Im (s™> M ' Sf > M ) ] P M (V)(7) 

Here G = e 1Tp t J l ps | JoVtip^sur is the elastic con- 
ductance of the junction while G$ = x 2 Gq, with 
x = Ts/Tq the relative intensity of the inelastic chan- 
nel. We define the spin polarization of the tip Vt = 
(PTT — p T jJ I (p T ^ + p T |) . Functions i± are defined as 
i±{A MtM >-eV) = Q(A MM >-eV)±g(A M ,M>+eV). The 
average adatom magnetization along the tip magnetiza- 
tion axis is (S z ) = Y^ai p m{V){M\S z \M) and depends 
on the bias voltage, the central result of the manuscript. 
Importantly, the magnetoresistive contribution to the 
elastic current permits to track changes in magnetiza- 
tion. Although x probably varies from system to system, 
we take x = 1, close to the results reported for a single 
Mn adatom [B| 

The inelastic contribution has also two terms, one of 
which depends on the tip polarization. The first inelastic 
term is proportional to i_, whose derivative with respect 
to V gives the characteristic inelastic conductance steps 
as \eV\ goes across ±Am,m'- By contrast, the Vt de- 
pendent term of the inelastic conductance, proportional 
to i+, yields steps at the excitation energies of opposite 
sign as the polarity of the bias is reversed. Both the elas- 
tic and inelastic term proportional to Vt can produce a 
dl/dV which is not an even function of bias. 

We now can address the main question of this 
manuscript: How does a spin polarized tunneling cur- 
rent affect the magnetization of a single adatom?. We 
consider first the simplest possible situation, where the 
magnetic moment of the tip is parallel to the easy axis 
of the single atom. We choose a single Mn atom on 
a CU2N surface, which in the case of non-magnetic 
tip is very well characterized experimentally [5| and 
theoretically d, [T3, EH, (3- The spin of the Mn atom 
in this situation is S = 5/2 with D = —0.039 meV and 
E = 0.007 meV. Since E « \D\ we can limit our qual- 
itative discussion to the case E = 0, so that the eigen- 
states of TLs are also eigenstates of S z . The numerical 
simulations are done with E ^ and do not change qual- 
itatively. In the absence of applied magnetic field and at 
temperatures much smaller than the zero field splitting 
4|Z?|, the equilibrium distribution is such that the two 
ground states, S z = ±5/2, are equally likely and the 
equilibrium average magnetization is zero. 

The non-equilibrium dynamics of the atom depends 
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FIG. 2: (a) Average magnetization (S z ) versus applied bias 
for T = IK (dashed line) and T = 0.1K (solid line), (b) 
Average magnetization after saturation (eV = 5meV) versus 
tip polarization for the same two temperatures, (c) dl/dV 
for single Mn in CU2N surface probed with a non-polarized 
tip (solid line) and two different tip polarizations, Vt = — 1/3 
(thin-dashed line) and Vt = — 1 (thick-dashed line) for a fixed 
temperature T — 0.1K. (d) Variation of the populations Pm 
with applied bias for non-polarized tip and T — 0.1K. 



on a number of parameters that can be tuned exper- 
imentally, such as the bias voltage, the temperature, 
and the ratio r = ^ LE -. The latter depends on the tip- 
adatom distance. As a general rule, the processes that 
drive the magnetic adatom out of equilibrium are pro- 
portional to Vti p wf UI whereas the processes that cool the 
spin down (provided that kbT < e\V\) are proportional 
to ftip+ v sur- Thus, the non-equilibrium effects are higher 
as r increases. We always take r < 1. Other parameters, 
such as the inelastic ratio x and the magnitude of the tip 
polarization Vt are not so easy to control. 

In Fig. El^a) we plot the steady state average magneti- 
zation of the Mn atom for Vt = — 1 and v s — 2vt = 1 
and two temperatures, T = IK and T = 0.1K, as ob- 
tained from solving Eq. ([3]). The tip polarization was 
assumed parallel to the Mn easy axis. The result has 
three outstanding features. First, the magnetization of 
the atom can be reversed from parallel to antiparallel just 
by electrical means [Fig. EJa)]- Second, the saturation 
magnetization can take very large values, which increase 
up to 100 percent for E = as the degree of spin po- 
larization of the tip increases [see Figl2fb)] . Third, the 
magnetization switch produces a large assimetry between 
positive and negative bias conductance G — dl/dV [Fig. 
[2fc)] which would be the experimental smoking gun of 
the spin transfer. 

The magnetization shown in Fig. 2(a,b) arises from the 
competition between spin-transfer from the spin polar- 
ized current to the atom and spin relaxation. The micro- 
scopic spin-transfer events, depicted in Fig. l(b,c,d) , are 
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the same than those resulting in steps in IETS [B|, 0| . In 
the case of Mn atom probed by non-magnetic tip, a single 
step has been reported [f| , related to the spin increasing 
transitions —5/2 — > —3/2 and its time reversal counter- 
part, the spin decreasing transition +5/2 — > +3/2, both 
with energy 4|D| (see Fig. ((Ifc,d)). For non magnetic 
electrodes there is no net spin transfer because the rates 
of these two processes are identical. In contrast, net 
spin transfer occurs when the tip is magnetic and these 
two processes have different rates. This is more easily 
seen in a extreme case Vt = 1, a half metallic tip such 
that pti — 0. Then, spin flip assisted tunneling from tip 
to surface can only decrease the transport spin (upper 
diagram in Fig. lb), increasing the spin of the magnetic 
adatom (process (I) and (II) in Fig. (lc)) whereas pro- 
cesses (III) and (IV) (Fig. lc) are not allowed, or for 
non-fully polarized tips, partly suppressed. When the 
bias is reversed, electrons flow from the surface to the 
tip. In this case, pt\ = implies that surface-tip spin- 
flip assisted tunneling can only increase the transport 
spin (lower diagram in Fig. lb), decreasing the adatom 
spin (processes (I) and (II) in Fig. Id). 

In the case of non-magnetic tip, our theory predicts a 
non-equilibrium effect which has been already observed 
experimentally jf|: the small decay of dl/dV for eV 
larger than the inelastic threshold. The spin-flip assisted 
events (Fig. l(b,c,d)) deplete the ground state doublet 
in favor of the first excited state doublet [fig. (J2Ji)]. This 
opens two new transport channels (±3/2 — > ±1/2, (pro- 
cesses (II) and (IV) in Fig. l(c,d)) which, in this case, 
happens to be less efficient than the channels ±5/2 — > 
±3/2 (processes (I) and (III) in Fig. l(c,d)). As the 
bias increases above 4|D|, the efficient inelastic transi- 
tions (I,III) are replaced by the less efficient transitions 
(II, IV), decreasing of the inelastic contribution to dl/dV. 

The spin transfer events that "heat" the adatom, re- 
quire that a quasiparticle has an excess energy larger 
than 4|D|. They compete with spin cooling events (order 
^tipi^sur) that are always available. At T = 0, the cur- 
rent induces changes in the magnetization only for bias 
larger than the inelastic threshold 4|£)|. At finite temper- 
ature there are always quasiparticles thermally activated 
above the 4|D| threshold and they are enough to induce 
the magnetization in the atom, even for eV < 4|_D|, pro- 
vided that a sufficiently long time is used in the process. 
In Fig. [3] we show how the time-scale necessary for the 
switching depends dramatically both on eV and k},T: the 
switching time can decrease by up to four orders of mag- 
nitude as either the bias goes above 4|£)| or the temper- 
ature is raised above 4|D|/fcf,. 

Finally, we have studied the effect of the deviations of 
the magnetization axis n with respect to the adatom easy 
axis. Curves similar to Fig. 2a are obtained for the av- 
erage magnetization of the adatom along ft as a function 
of eV, but the saturation magnetization at high bias is 
a decreasing function of the angle formed by n and the 



adatom easy axis. Whereas a small tilt (5 degrees) barely 
changes the curve compared to Fig. 2a, the efficiency of 
the spin transfer is minimal when ft is perpendicular to 
the adatom easy axis but still can reach 0.6. 
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FIG. 3: Switching time TVei for Vt = —1, vs = 1 and vt = 
0.5. (a) T re i as a function of applied bias V for T = 0.1K and 
(b), as a function of temperature for eV = 2\D\. 

In conclusion we have shown that the spin of a single 
magnetic adatom can be polarized, reversed and moni- 
tored by means of spin polarized currents in the STM 
configuration. The time scale for the single spin switch- 
ing can be as quick as a few nanoseconds. Our pro- 
posal adds the possibility of single spin control to the 
wide range of uses of STM in the field of nanospintron- 
ics 0, This possibility is based on non-equilibrium 
processes and differs from recent work in which switch- 
ing is achieved by a change of sign of the equilibrium 
tip-adatom exchange interaction [191 ]. 
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